Draft version November 11, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



TRACERS OF CHROMOSPHERIC STRUCTURE I: OBSERVATIONS OF CA II K AND Ha IN M DWARFS 

Lucianne M. Walkowicz 

Astronomy Department, University of Washington, Box 351580, Seattle, WA 98195-1580 



00 

o 

O 

> 
O 



43 

6 



> 

oo 
o 



oo 
o 



X 



Suzanne L. Hawley 

Astronomy Department, University of Washington, Box 351580, Seattle, WA 98195-1580 
Draft version November 11, 2008 

ABSTRACT 

We report on our observing program^] to capture simultaneous spectra of Ca II and Balmer lines in 
a sample of nearby M3 dwarfs. Our goal is to investigate the chromospheric temperature structure 
required to produce these lines at the observed levels. We find a strong positive correlation between 
instantaneous measurements of Ca II K and the Balmer lines in active stars, although these lines may 
not be positively correlated in time-resolved measurements. The relationship between Ha and Ca II K 
remains ambiguous for weak and intermediate activity stars, with Ha absorption corresponding to a 
range of Ca II K emission. A similar relationship is also observed between Ca II K and the higher order 
Balmer lines. As our sample consists of a single spectral type, correlations between these important 
chromospheric tracers cannot be ascribed to continuum effects, as suggested by other authors. These 
data confirm prior non-simultaneous observations of the Ha line behavior with increasing activity, 
showing an initial increase in the Ha absorption with increasing Ca II K emission, prior to Ha filling 
in and eventually becoming a pure emission line in the most active stars. We also compare our optical 
measurements with archival UV and X-ray measurements, finding a positive correlation between the 
chromospheric and coronal emission for both high and intermediate activity stars. We compare our 
results with previous determinations of the active fraction of low mass stars, and discuss them in 
the context of surface inhomogeneity. Lastly, we discuss the application of these data as empirical 
constraints on new static models of quiescent M dwarf atmospheres. 

Subject headings: stars: late-type — stars: low-mass, brown dwarfs — methods: data analysis 



1. INTRODUCTION 

Radiation from the chromospheres of Sun-like stars 
comes primarily from the Ca II resonance lines in the 
blue and Ca II infrared triplet, the Fe II and Mg II res- 
onance lines in the nea r UV, and the Balmer series of 
Hydrogen in the optic al (jLinskv et al. 1982; Fouka1 [l99d : 
iReid fe; Haw ley 200^). In M dwarfs, high surface grav- 
ity compresses the stellar atmosphere, leading to higher 
density c hromospheres and particularly strong Balmer 
emission (Lin skv et al.lll982ft . While the Ca II lines are 
the brightest emission lines formed in the of the Solar 
chromosphere, contributing a radiative loss of ~ 5xl0 6 
ergs cm" - 2 s- 1 (|Foukallll99a ). the Balmer emission in M 
dwarfs is much larger than all other chromospheric lines, 
compri sing ^10~ 4 of the bolometric lumino sity at early 
types (jLinskv et alJll982t IWest et alJ 12004ft . This fact, 
combined with the intrinsic faintness of M dwarfs in the 
blue, has made the red Ha line the primary diagnostic 
spectral feature of chromospheric activity in these stars. 
Many M dwarfs also have prominent Ca II resonance line 
emission- particularly Ca II K, which is unblended even 
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at low resolution (as opposed to Ca II H, which is blended 
with He). The subset of these stars which also show Ha 
emission are designated "active" (dMe). 

Ha and Ca II K are two of the strongest optical emis- 
sion lines in active M dwarf chromospheres, responsible 
for cooling the atmosphere, balancing the magnetic heat- 
ing, and determining the resulting equilibriu m structure 
(jSundland et alj!988tlMauas fe Falchilll994f ). Across the 
M spectral class, there is a range of emission strength in 
Ca II K, and a wide variety of both absorption and emis- 
sion in Ha. Ha appears to trace hotter regions of the 
chromosphere (> 7000 K), while Ca II K traces the cooler 
reg ions between the temperature minimum and ~ 6000 
K (iGiampapa. Worden. fe LinskvH l982t IC ram & Mullan 
ll98a IWalkowi cz et al.ll2009ft . Thus together, Ha and Ca 
II K offer complementary information on chromospheric 
structure. 

The nature and relationship of the Ca II K and Ha 
lines in both active and inactive M dwarfs has been con- 
troversial. It may be that Ca II K and Ha are not only 
complementary but integral pieces of information, which 
separately cannot offer a true understanding of the pro- 
cesses they trace. One point of contention has been the 
"zero point" of chromospheric activ ity, or what defines 
the weakest existing chromospheres. lYoung et al.l (|1984ft 
have argued that the zero point for chromospheric ac- 
tivity is defined by strong H a absorption. In contrast, 
Stauf feTfe Hartmannl (|1986ft claim that those stars with 
the weakest chromospheres also have the weakest Ha ab- 
sorption. In their scenario, Ha begins as a weak, purely 
photospheric absorption line, which (when the atmo- 
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sphere is subject to magnetic heating) first increases in 
absorption equivalent width, then decreases as it fills in 
with emission, then becomes a pure emission line. The 
implication of this scenario is that both stars with in- 
trinsically weak chromospheres, and stars with moder- 
ate amounts of magnetic activity, may exhibit weak Ha 
absorption. Without additional chromospheric indica- 
tors, such as the Ca II K line, it is impossible to dis- 
tinguish intermediate activity stars from inactive stars. 
In M dwarfs, the cold tem perature minimum is largel y 
transparent to Ha photons (|Pettersen fe Colemanlll98lf ). 
so only a very w eak Ha absorption lin e is produced in 
the photosphere (|Cram fc Mullanl[l985h - thus we adopt 
those stars having weak Ha absorption (< O.lA) and 
weak Ca II K emission (< O.lA) as being representative 
of the "zero point" of activity in our sample. 

Prio r studies of the relationship between Ha and Ca 
II K (IRobinson et all 119901: iRauscher fc Marcvi 120061 : 
ICincunegui et alJ 120071 ) have routinely observed a pos- 
itive correlation between these lines i n activ e M 
dwarfs. Howev e r bo th IRobinson et alJ (|1990l ) and 
ICincunegui et al.1 (|2007| ) have suggested that this ob- 
served correlation is due to a color effect- in other words, 
that the relationship between these lines observed in a 
sample of mixed spectral types is due to enhanced con- 
tribution from high continuum fluxes in earlier type stars. 
These authors argue that because time-resolved observa- 
tions of Ca II K and Ha emission are not always posi- 
tively correlated in individual stars, the positive correla- 
tion between instantaneous measurements of these lines 
is the result of their spectral type, rather than the prop- 
erties of their chromospheres. 

Observations of Ca II K and Ha provide an impor- 
tant "point of contact" in linking theory to actual M 
dwarfs. Previous efforts to study this relationship have 
often relied on observations from different epochs to 
provide measuremen ts of both lines for a give n star 
ijCram fc Giampapa|[l987l : iHoudebine fc Stempelslll997f ) . 
In active M dwarfs, however, both Ca II K and Ha vary 
on timescales from minutes to decades. Therefore, mea- 
surements of lines obtained from non— simultaneous ob- 
servations are not truly comparable, and any attempt 
to model them may shed no greater light than compar- 
ing lines in two totally different stars. Indeed nearly 
every previous study, observational and theoretical, cites 
non— simultaneous data as the chief source of uncertainty 
in observed relations, and reco mmends synoptic observa- 
tions of a sample of M dwarfs (|Cram fc Giampapalll987l : 
IRobinson et~atlll990f) . 

Although many surveys have measured the Balmcr 
lines of low mass stars, particularly Ha (e.g. the 
Palomar-Michigan State Un iversity Survey, or PMSU; 
iGizis. Reid. fc Hawlevll2002t ). the desired dataset of syn- 
optic Balmer and Ca II K observations is sparse. The 
non-simultaneous observations taken in previous surveys, 
along with samples that span a wide range of M spectral 
types (and therefore a wide range in color and absolute 
magnitude), have left the relationship between these im- 
portant chromospheric tracers ambiguous. To remedy 
this uncertainty, we have carried out an extensive observ- 
ing program of Ha and Ca II K in nearby M3 dwarfs. By 
taking simultaneous data for a single spectral class, we 
are able to determine if the observed correlation between 



these chromospheric tracers is due to a continuum effect, 
or to an actual relationship between these lines in active 
chromospheres. We also examine whether multiple ob- 
servations of the Balmer and Ca II lines in an individual 
star are correlated over time in the same sense as single 
observations, and if not, how this behavior can be recon- 
ciled with instantaneous measurements of these lines in 
single observations of that same star. 

We also compare these observations with previous de- 
terminations of the act ive fraction of low mass stars 
(|West et al.ll2004l2008al ). which have relied solely on the 
Ha line as an activity proxy, although it has been widely 
assumed that most M dwarfs possess some level of Ca 
II K emission. Using the Ca II K line as an additional 
activity measure, this survey will allow us to determine 
what fraction of stars deemed "inactive" on the basis of 
Ha absorption actually do have Ca II K emission and 
therefore belong to an intermediate activity population. 

Additio nally, we examine th e prop osed Ha behavior 
found by lStauffer fc Hartmannl (|1986l ). whereby Ha first 
increases in absorptio n depth prior to filling in an d go- 
ing into emission. The IStauffer fc Hartmannl (|1986f ) sce- 
nario was derived from non-simultaneous observations 
of a sample com posed of cool stars of a relatively wide 
range of mass. ICram fc Mullanl (|1985h reported some 
success reproducing these observations using very simple 
model chromospheres, consisting of a slab with a con- 
stant temperatu r e and density. Using this simple model, 
I Cram fc Mullanl f| 19851 ) predict a different slope for the 
filling in of the Ha line depending on whether the in- 
crease in equivalent width is due to an increased heating 
rate in a homogeneous chromosphere, or an increased fill- 
ing factor of active regions in a laterally inhomogeneous 
chromosphere. If the filling in of the Ha line is evident 
in our observations, it will provide further constraints on 
how chromospheric structure changes with increased ac- 
tivity, as well as the role played by lateral inhomogeneity. 

We are motivated to investigate s urface inhomogene- 
ity by the results of We st et all (|2008aD . who used 
the observed fraction of active stars per spectral type, 
along with a simple dynamical model of the Galaxy, 
to infer activity lifetimes for M dwarfs across the spec- 
tral class. A large increase in the activity lifetime 
occurs over the range of M3-M5, corresponding to 
the mass boundary where stars are believed to be- 
come fully convective (~ 0.35M W ; iReid fc Hawlevl 120051: 
IChabrier fc Baraffdfl997h . The transition to full convec- 
tion may mark a change from a Solar-type, rotationally- 
dependent magnetic dynamo, to a turbulent dynamo 
that may or may not depend strongly on rotation. The 
decline in activity in a given star over its lifetime may 
also indicate a change in how that star's activity is gener- 
ated: objects of transitional mass (such as the M3 stars 
we examine here) may still possess some marginal radia- 
tive zone, though the majority of the star is characterized 
by deep convection. In this case, it is possible that during 
the star's youth the activity is dominated by a Solar-type 
dynamo, where the boundary of the convective and radia- 
tive zones is responsible for most of the field generation, 
though a turbulent dynamo is also at work. As the star 
spins down with age, the rotation-dependent aQ dynamo 
declines, and the turbulent dynamo becomes the predom- 
inant means of field generation. Indeed, recent simula- 
tions of the magnetic dynamo in fully convective stars 
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by iBrowniiiel (|2008f ) show that the turbulence-generated 
field has both large- and small-scale features. The typi- 
cal size of these field structures are related to the radius, 
with those generated deeper in the interior being larger 
than those near the surface. In addition, the large scale 
field in these fully convective models can be generated 
and sustained even though differential rotation (which 
plays such a key role in Solar-type dynamos) is eventu- 
ally quenched by the field itself. Interestingly, although 
the iBrownind (|2008l ) models are convective throughout, 
convection in the core is relatively weak, and most of the 
energy is transported via the radiation. 

There is some precedence for this "dual-dynamo" sce- 
nario on the Sun, which possesses both a large-scale 
global field and a small scale magnetic component that 
seems to be independent of the global field. The gross 
magnetic features of the Sun are thought to be due 
to the afl dynamo that originates in the tachocline, 
but the magnetic field of the quiet Sun is likely gener- 
ated by tur bulent convective motions close to the stel- 
lar surface (iMeneguzzi fc Pouauetl Il986t iDurnev et al.1 
119931 I Voder fc Schusslerl |2007ft . The chromosphere is 
a thin layer, dynamic but not convective, threaded by 
convection-twisted mag netic fields th at form the quiet 
Sun magnetic network (Abbett 2008). Field generation 
by turbulence near the surface is not to be confused with 
theories of a 2 dynamos, which generate global scale tur- 
bulent fields, or with acoustic heating of the chromo- 
sphere. Although acoustic heating is bound to be at 
work wherever convective motion is involved, purely me- 
chanical means apparently cannot fully account for the 
chromospheric and coronal emission in the Sun, and the 
heating due to these surface turbulent fields is most likely 
magnetic in nature (jBercik et al.ll2005l . 1200(1 . 

In the stellar atmosphere, the manifestation of this be- 
havior may be surface inhomogeneity, where the chro- 
mosphere can be represented by two components: one 
corresponding to magnetically active plage regions re- 
lated to the large-scale field, analogous to active re- 
gions on the Sun, the other component representing a 
basal , turbulence-driven chromosphere (Sc hrijver et al.l 
1989). Observations of the solar surface indicate that 
these inhomogeneities are complex indeed- Ca II K 
core emission corresponds spatially to regions of con- 
centrated magnetic field, such as in active plage regions 
and bright network grains, while Ha absorption is pro- 
duced in the both the mass of fibrils protruding from 
active regions, in mottles across the network of the quiet 
Sun and i n enhanced emission from bright points dur- 
ing fl ares (|Hasan fc van Ballegooijenll2008l : lRuttenl[2006l 
2007). The diversity of structures in the solar chromo- 
sphere seems to indicate that even a dual component 
description is somewhat simplistic in nature- the mag- 
netic features of a given star likely form a more contin- 
uous distribution, where t he distinction b etween active 
and quiet areas is unclear (|Schri ivcr 1988). In the stel- 
lar case, however, we have no such detailed observations 
of surface inhomogeneities and must proceed cautiously 
with the number of free parameters we allow in our in- 
terpretation of spatially averaged spectra. However, if 
the Ca II or Balmer lines preferentially trace particular 
chromospheric components, examining their relationship 
may shed light on the nature of magnetic structures and 
the stellar dynamo in convective-boundary mass stars. 



In the next section, we describe the observed sample. In 
Section 3, we discuss the observations and analysis, with 
particular attention to the effect of spectral resolution on 
interpreting relationships between chromospheric lines. 
In Section 4, we present the observed relationships be- 
tween Ca II and the Balmer lines. We also place these ob- 
servations in the context of previous surveys, addressing 
their implications for chromospheric structure. Lastly, 
we comment on our proposed theory of atmospheric in- 
homogeneity in light of the observations, and discuss fu- 
ture plans to test this proposition with a new suite of 
non-LTE model atmospheres. 

2. THE SAMPLE 

Our sample consists of 81 stars in the Palo- 
mar/MSU Nearby Star Spectroscopic Survey (PMSU ; 
Reid. Hawlev fc Gizis|[l995: Hawlev, Gizis & Reid] 11996 
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of the same approximate spectral type, M3V (see Table 1 
for a summary of our sample stars). The PMSU Survey 
determined spectral types b ased on a set of spectroscopic 
indices for TiO and CaH fsee lReid. Hawlev fc Gizill995l 
for details). Our sample is also quite homogeneous in 
color, having a median Cousins V-I = 2.55 ± 0.12 (cor- 
responding to a differ ence in T e ff of roughly ± 50K; 
iReid fc Hawlevl pOOl ). The PMSU Survey Ha mea- 
surements enabled us to select a sample that is represen- 
tative of varying levels of chromospheric activity, from 
very weak Ha absorption (indicating either a weak or 
intermediate strength chromosphere), to powerful emis- 
sion (indicating strong chromospheric activity). Limit- 
ing the spectral ran ge of our objects is important, as it 
has b een been noted (|Stauffer fc Hartmannlll986l : iRuttenl 
1987) that the relation between the chromospheric flux in 
the Ca II K line and the equivalent width in Ha depends 
on the spectral type of the star (i.e. the stellar contin- 
uum flux). Previous studies have attempted to correct 
for th is effect by binning t heir measurements by stellar 
color ([Robinson et al.l fl990L By limiting our sample to a 
single spectral class, we are able to examine the relation 
between Ca II K and Ha without having to disentangle 
systematic trends with color and spectral type. Although 
many of the stars we observe have existing optical data, 
our new data comprise simultaneous observations of both 
lines at comparable resolution. 

3. OBSERVATIONS AND ANALYSIS 

Data were obtained with the ARC 3.5m Telescope 
at Apache Point Observatory (APO), using the Dual 
Imaging Spectrograph (DIS) and the ARC echelle spec- 
trograph (ARCES). Additional high resolution data 
(HIRES) were available for many of our stars from 
the California Pla net Search (CPS; iButler et all 119961 : 
IWright et aLll2004D . While the DIS spectra provide reli- 
able flux information for our lines of interest, the echelle 
spectra enable an examination of the detailed line pro- 
file. High resolution data are particularly important for 
identifying intermediate activity objects, whose Ca II K 
chromospheric emission cores may not be resolved in the 
lower resolution spectra. Figure [1] shows example HIRES 
and DIS spectra for three stars of different activity lev- 
els. In the highly active star (top row: AD Leo), the 
Ha and Ca II K lines are quite clearly in emission in 
both the echelle and low resolution spectra, indicating 
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the presence of substantial chromospheric heating. The 
moderate activity example (middle row: Gl 362) shows 
very little Ha emission, but still has a noticeable Ca II 
K emission line. In the weakly active example (bottom 
row: GJ 1125), the Ha line is in absorption and Ca II K 
emission is not evident in the low resolution spectrum. 
Based on its Ha line, this star would be classified as "in- 
active" . However, the echelle spectrum of this star illu- 
minates the importance of high resolution observations- 
even the "inactive" object possesses a noticeable Ca II 
K emission core when examined in detail. In the fol- 
lowing subsections, we describe the observations, data 
reduction, and method of line measurement for the DIS 
and echelle spectra. We also compare high and low res- 
olution measurements of the same stars, to investigate 
what biases may be introduced by the choice of spectral 
resolution. 

3.1. DIS spectra 

DIS is a low dispersion spectrograph with two cameras 
fed by a beam splitter, allowing blue and red spectra to 
be obtained simultaneously. Our data were taken using 
a l'/5 slit in combination with the higher resolution grat- 
ings, which have a dispersion of 830.8 lines mm -1 on the 
red side (giving a dispersion of 0.84A pixel -1 , R = A/AA 
~ 3200) and 1200 lines mm -1 on the blue side (giving 
a dispersion of 0.62A pixel -1 , R ~ 2000). Spectra typi- 
cally covered the wavelength regions of 3700A to 4400A 
in the blue and 5600A to 7100A in the red. 

The DIS data wer e reduced using stand ard IRAF re- 
duction procedures (jMassev et al. I Il992t ). Data were 
bias-subtracted and flat-fielded, and one-dimensional 
spectra were extracted using the standard aperture ex- 
traction method. A wavelength scale was determined 
for each target spectrum using HeNeAr arc lamp cali- 
brations. Flux standard stars from lOkel (|1990l ) were ob- 
served each night and used to place the DIS spectra on 
a calibrated flux scale. 

3.2. Echelle Spectra 

The ARC echelle spectrograph (ARCES) captures the 
entire visible spectrum (roughly 3800A to 9800A) with 
a resolution of R~31,500. While this resolution is high 
enough that even small Ca II K emission cores can be 
measured from our spectra, it does not fully resolve the 
Ca II K line profile, and thus information regarding the 
presence or lack of a central reversal in the detailed line 
profile (an important tracer of NLTE effects in the upper 
chromosphere) is lost in these spectra. Cosmic rays were 
removed from the raw echelle images, which were then 
bias-subtracted and flat-fielded. Scattered light in the 
ARCES spectra was removed to the 6-8% level using the 
apscatter task in the IRAF echelle package, which fits 
a two-dimensional surface to the dispersed and cross- 
dispersed interaperture light. One-dimensional spectra 
were extracted for data and ThAr arc lamp images. The 
ThAr arc lamp spectra were then used to determine the 
wavelength scale for each target spectrum using ThAr 
arc lamp spectra. 

The CPS spectra were taken using the HIRES spec- 
trograph on Keck I (R-57,000; IVogt et al.lfl994h . and 
cover the wavelength regime of 3645-7990A. The Keck 
HIRES data were reduced using the Planet Search Re- 



duction Pipeline, describ ed in Wright et al. 2004 and 
iRauscher fc Marcvl (|2006[ L The pipeline removes scat- 
tered light from the raw spectra by fitting a smooth func- 
tion to counts lying between the orders, and subtract- 
ing this level from the raw data. This method removes 
the scattered light to ± 1 photon (G. Marcy, private 
communication). Equiv alent widths for these ob jects 
have also be en rep orted in lRauscher &; Marcvl (|2006l ) and 
IWest et al.1 (|2006h . but we have remeasured them as part 
of our analysis to ensure consistency. 

3.3. Measurement of Ca II K and Balraer Lines 

The observed line profiles are a superposition of the 
photospheric and chromospheric contributions. In the 
case of the Ca II K line, the chromospheric emission core 
is easily identifiable in our echelle spectra even for stars 
with low activity levels, and so the two components of the 
line can be disentangled. As the distinction between pho- 
tospheric and chromospheric contribution to the Balmer 
lines is considerably more ambiguous, we measure the en- 
tire (photospheric+chromospheric) line. Measuring the 
entire Ha line enables direct comparison between our 
work and previous studies that use Ha as an activity in- 
dicator, where traditionally no distinction has been made 
between the chromospheric and photospheric contribu- 
tions to the line. 

The Ca II K line was measured in the same way 
for both the HIRES and ARCES echelle data as well 
as the low resolution DIS data. As the Ca II K line 
profile contains both contributions from the chromo- 
spheric emission and the photospheric absorption, it 
is necessary to correct for the photospheric contribu- 
tion in order to study the chromospheric emission it- 
self. In order to remove the photospheric contribution 
to each line, we use a radiative equilibrium NextGen 
model atmosphere (log g=5, T=3400K, solar metallicity; 
lHauschildt et all 119991) with the NLTE radiative trans- 
fer code 'RH' (jUitenbro ek 2001) to calculate the photo- 
spheric Ca II K line profile. The calculated line profile 
was normalized and scaled to match the continuum near 
the Ca II K line for each star. The model photospheric 
line profile was then overplotted with the observed data, 
and a line integration region was interactively chosen to 
include only the chromospheric emission; this method is 
illustrated in Figure 4. As the stars in our sample are all 
the same spectral type, the same model profile is used in 
each case and the correction is applied in the same way 
to each spectrum. 

The total area under the interactively chosen region is 
then integrated both for the observed spectrum and the 
radiative equilibrium model. The flux above the radia- 
tive equilibrium model is taken as the chromospheric flux, 
and divided by the mean continuum level measured from 
a relatively featureless portion of the nearby spectrum to 
yield the equivalent width of the line (a summary of the 
wavelength regions chosen for the mean continuum and 
line integration regions is given in Table 2). We com- 
pare the results obtained with photospheric correction 
for objects with both echelle and DIS spectra in the next 
section. 

As the Balmer lines have relatively simple profiles com- 
pared to the Ca II K line, such a detailed method of mea- 
surement was not necessary. They were instead measured 
using an automated method that sums the flux in rectan- 
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gular passbands, using a fixed integration region for each 
line. The Balmer lines in the DIS spectra were measured 
using the same fixed integration and continuum regions 
as for the high resolution observations (specified in Table 
2). 

As the lower resolution DIS spectra are flux calibrated, 
fluxes for the lines could be measured in addition to 
equivalent widths. Where given, surface fluxes were cal- 
culated using distances from PMSU (derived from the 
Mt,-Ti05 spectroscopic parallax relation, and calibrated 
against those stars in the sample posessing trigonometric 
parallaxes; see lGizis. Reid. fc Haw lev 200 21) and a typical 
radiu s for an M3V star (2.85 xlO 10 cm: iReid fc Hawlevl 
I20051) . 

Table 3 provides the line fluxes and equivalent widths 
for the entire sample as measured from the low resolution 
DIS spectra, while Table 4 provides equivalent widths 
measured for the subset of stars with echelle observa- 
tions. 

3.4. Comparison Between Low Resolution and Echelle 

Data 

The detailed profiles from our echelle spectra enable us 
to easily disentangle the photospheric and chromospheric 
components of our lines, which is particularly important 
for stars with low levels of activity (and thus ambiguous 
low resolution spectra). We use the comparison between 
the equivalents widths measured from the low and high 
resolution spectra with photospheric correction to cali- 
brate the effect of spectral resolution on our results. 

Figures [2] and [3] compare the equivalent widths ob- 
tained from the DIS spectra with those from the echelle 
spectra. For low levels of activity, the Ca II K equivalent 
widths derived from the corrected low resolution data 
are a decent match with those derived from the high res- 
olution data. The difference in equivalent widths calcu- 
lated from the ARCES spectra versus the DIS spectra are 
±0.70A in Ha and ±l.llA in Ca II K, while equivalent 
widths measured from the HIRES spectra are ±0.17A 111 
Ha and ±0.18A in Ca II K. However, as the activity of 
the star increases, equivalent widths calculated from the 
low resolution data are a good deal higher in both Ca II 
K and Ha: in Ha, the equivalent widths measured from 
ARCES spectra of active stars are 2.13A lower than those 
measured from the low resolution spectra, while HIRES 
measurements are 0.99A lower. In Ca II K the effect is 
more pronounced, with the equivalent widths measured 
from ARCES and HIRES spectra being 1.63A and 3.73A 
lower, respectively. 

As our echelle spectra are not flux calibrated, the con- 
tinua have subtle differences in slope compared to the 
continua of the DIS spectra. While this effect is small for 
the majority of the sample, which have small Ha equiv- 
alent widths, it becomes much more noticeable when the 
line is in emission and contains considerably more flux. 
In the case of the HIRES spectra, Ha falls very close to 
the edge of the reddest order, and there is no adjacent or- 
der with which to join it. In both the HIRES and ARCES 
spectra, the equivalent widths calculated from the echelle 
spectra are smaller than those calculated from the DIS 
spectra. 

The differences between the two datasets indicate that 
equivalent widths calculated from high resolution data 
will be smaller for both the Ha and Ca II K lines in ac- 



tive stars. Therefore, measurements taken from the high 
resolution spectra will be offset in their absolute values 
of the equivalent widths for the same star. However, 
comparisons between the high resolution measurements 
of Ha and Ca II K will show the same relationship as 
for the low resolution data due to the net effect being 
in the same direction for both lines. If correction for 
the photospheric component of the Ca II K line is not 
performed, however, the activity level will always be un- 
derestimated, especially for low activity objects whose 
emission cores are not resolved. In the following section, 
we only compare data of comparable resolution. 

4. RESULTS 

In this section we present our results and interpret 
them in response to the questions raised in the Intro- 
duction. 

4.1. The Relationship Between Ha and Ca II K with 
Increasing Activity 

Comparison of the Balmer and Ca II line s reve als a 
number of interesting relationships. Figure 5(a) plots 
the equivalent width of Ha versus the equivalent width 
of Ca II K as measured from our DIS spectra. The Ca 
II K line has been corrected for the photospheric compo- 
nent as described in Section 3, so the equivalent widths 
cited here represent the purely chromospheric contribu- 
tion to the line. For the most active stars (those with Ha 
in emission, indicating strong chromospheric heating), 
the Ha and Ca II K equivalent widths arc clearly corre- 
lated and increase monotonically, with increased scatter 
for stronger line strength. A linear fit to the Ha and 
Ca II K equivalent widths for the active stars yields an 
RMS deviation of ±0.38A, with stars of Ca II K equiv- 
alent width of less than ^5.35A X having an RMS devi- 
ation of ±0.2lA, while the most active stars (EWcaiiK 
> 5.35A) have an RMS deviation of ±0.56A. As these 
data are measured from a homogeneous sample of stars 
of the same spectral type, and the observations of Ca 
II K and Ha are simultaneous, the increased scatter is 
likely intrinsic and suggests a mild decoupling of the line 
behavior in the strongest chromospheres. Stars possess- 
ing weak to intermediate chromospheres cluster in the 
lower left corner of the plot, and are spread over a wide 
range in Ca II K for a relatively small range in Ha. These 
wea kly a ctive stars, plotted at an expanded scale in Fig- 
ure |5(b)| do not show any obvious correlation between 



the Ha and Ca II K equivalent widths. That a range in 
Ca II K emission strength can be associated with cither 
deep or shallow Ha ab s orptio n lines lends credence to the 
Stau ffer~fc Hartmannl (|1986l ) scenario, where the photo- 
spheric Ha absorption line "fills in" with chromospheric 
emission as activity increases. However, an initial deep- 
eni ng of Ha absorption as C a II K increases, suggested 
bylCram fc Giampapal (|1987l ) and IStauffer &: Hartmannl 
(1986), is not obvious in these low resolution observa- 
tions. 

This value corresponds to the mean of the Ca II K equivalent 
widths measured for active stars and is a largely arbitrary choice for 
the division between "active" and "very active". Our result that 
the scatter increases with increasing activity is fairly insensitive to 
this value, routinely yielding an RMS deviation for the most active 
stars that is roughly twice that for less active stars. 
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Figure 5(b)| illustrates the number of stars that would 
traditionally be classified as "inactive" based on having 
Ha in absorption, which actually do possess some chro- 
mospheric activity (indicated by their Ca II K emission) . 
Of the 81 stars in our sample, 46 stars were classified 
as inactive based on having Ha in absorption; with the 
additional information provided by our Ca II K measure- 
ments, we find that all but possibly 1 of these stars have 
at least a low level of Ca II K emission. As discussed 
in the introduction, this result implies that Ha-based 
measurements of the active fraction in early M dwarfs 
(above the fully-convective mass boundary) may trace 
what fraction of stars have a solar-type magnetic compo- 
nent that relics on the presence of a radiative-convective 
zone boundary (the tachocline). However, as Ha emis- 
sion is not only present but increasingly common in very 
low mass, fully convective M dwarfs which do not possess 
a tachocline, Ha in these stars may be associated with 
a global turbulent dynamo and thus trace an altogether 
different source of magnetic activity (as discussed in the 
Introduction) . 



Figure 6(a) shows echelle equivalent widths of the Ca 
II K chromospheric emission cores against the full (pho- 
tospheric+chromospheric) equiv alent width o f the Ha 
line as a comparison with Figure 5(a) Figure 6(b) pro- 
vides a simi lar comparison to the low ac tivity stars as 



Figure 5(b) 



Figure 6(a) shows a 



As in Figure 5(a) 
very small variation in H a abs orption for a range of 
Ca II K emission. Figure 6(b) is intriguing, however, 



as the U-shape locus of points does seem to indicate 
an initial deepening of the Ha absorption line prior to 
its filling in and going i nto emission. This relation- 
ship was also observed by iRauscher fc Marcvl (|2006f ) in 
their analysis of Ca II K HIR ES data and Ha equiv- 
alent widths from iGizis. Reid. fe Hawlevl (|2002t ) for M 
dwarfs of mixed spectral types, but the trend persists 
if we examine a homogeneous sample. Arrows indicat- 
ing the approximate place ment of the locus predicted by 
ICram fc Giampapal (|1987l ) are shown with the data to 
guide the eye. The maximum Ha absorption equivalent 
width is ~-0.35A, which is consi s tent w ith the range ob- 
served by [StaufeTfeHirtmin^ (HHl) . The maximum 
in Ha absorption seems to occur at roughly 0.4 - 0.5A Ca 
II K equivalent width, suggesting that stars having Ca 
II K emission at this level and above have passed some 
critical threshold for the onset of Ha emission. 

4.2. Chromospheric Coupling: Surface Fluxes and 
Higher Order Balmer Lines 

In Figure [3 the Ha equivalent width is plotted against 
the chromospheric Ca II K surface flux for our DIS sam- 
ple. For M stars with Ha in absorption, a relatively 
small variation in Ha equivalent width corresponds to a 
wide ran ge of Ca II K emission, confirming the results 
of Figure 5(a) above. For the strongest chromospheres, 
the de coupli ng of Ca II K and Ha emission suggested by 
Figure 5(a) is much more evident, with strong Ca II K 
flux being associated with a wide range of Ha emission 
streng ths. These results confirm those of iRobinson et all 
(1990), although their analysis was complicated by the 
spread in spectral type of their sample (early K to mid 
M). They similarly corrected the Ca II K for the photo- 
spheric contribution to the equivalent width, but used a 
solar model scaled to the various temperatures of their 



sample stars. 

Figure [8] shows the relationship between the equiva- 
lent width of Ca II K a nd tho se of H<5, H7, He and H8. 
As in the case of Figure 5(a) the equivalent widths are 
calculated from single DIS observations of our sample 
stars. Our result s provide an interesting comparison to 
IWest et all (|2006| ). who reported that in multiple obser- 
vations of the Balmer and Ca II lines in a given star, 
the higher order Balmer lines and Ca II H and K lines 
trace one another, while Ha varies seemingly randomly. 
From this plot, it is evident that the higher order Balmer 
lines follow Ca II K in a simil ar way to Ha. In co ntrast, 
Figure M (reprinted here from IWest et ail l2008bh shows 
the Balmer and Ca II K line equivalent widths measured 
from repeat echelle ob servations o f AD L eo. Confirm- 
ing the observation of We st et aTl (|2006f ). we see that 
while the higher order Balmer lines vary little in rela- 
tion to Ca II K (AEW H(5 =±0.16A; AEW ff7 =±0.15A; 
AEW ffe =±0.12A), the Ha equivalent width can vary 
dramatically o ver ti me (AEW# Q =±0.57A). Taken to- 
gether, Figures [5(a)| El and [9] imply that while the emis- 
sion in Ca II K and the Balmer lines are positively corre- 
lated when measured from a single observation of a given 
star, they are not necessarily positively correlated over 
multiple observations of that star. In other words, more 
active chromospheres produce greater emission in both 
Ca II K and the Balmer lines overall, but Ha and Ca 
II K do not always trace one another directly in time- 
resolved measurements. 

4.3. Investigation of Continuum Effects 

Although our sample is comprised of stars of the same 
spectral type, our stars cover a range of ~2 magn itudes 
in M„ . As previous stud ies bv lRobinson et all (jl990f ) and 
ICincunegui et all (|2007f ) have suggested that the positive 
correlation between Ca II and the Balmer lines in active 
stars is due to enhanced continua in early spectral types, 
we investigated the possible relationship between abso- 
lute magnitude and activity indicators for our sample. 
Figure [10] shows the continuum fluxes near Ha and Ca 
II K versus the absolute V magnitude, while Figure QT] 
shows the equivalent widths of Ha and Ca II K versus 
absolute V magnitude. The absence of any relationship 
between absolute magnitude and the level of either activ- 
ity or continuum flux shows that the positive correlation 
we observe between Ha and Ca II K in active stars can- 
not be ascribed to a continuum effect. 

4.4. Interpretation 

A possible interpretation of these results is that the Ca 
II K and Ha fluxes observed are due to a laterally inho- 
mogeneous atmosphere. Stated simply, increased activ- 
ity leads to stronger emission in all lines, but does not 
imply that these lines trace the same source of activity 
or the same regions o n the star. In analogy t o the "zebra 
effect" suggested by iPettersen et al.l (|1992T ) to describe 
photometric variation on BY Dra stars, one may picture 
the stellar atmosphere as consisting mostly of a cooler 
'basal' chromosphere, heated either acoustically, through 
a weak, pervasive magnetic field, or some combination of 
the two. Stars with greater activity may have, in ad- 
dition to this basal chromosphere, hotter active regions 
where magnetic flux tubes thread through the stellar at- 
mosphere, analogous to bright plage regions on the Sun. 
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The stronger the stellar magnetic field, the more area 
of the stellar surface is covered by these active regions. 
The observed line profiles are a superposition of the line 
profiles produced by these two atmospheric components. 

The Ca II K emission core responds linearly with 
increased activity- more active chromospheres produce 
stronger Ca II K emission. The Ha line, on the other 
hand, may be found in various degrees of absorption for 
low to intermediate levels of activity, and in emission for 
strong activity. Let us describe a few limiting cases: (1) 
the majority of the stellar surface is covered with quiet 
chromosphere and a small fraction of very active regions; 
(2) the majority of the stellar atmosphere consists of the 
active component, with a small fraction of the quiet chro- 
mosphere, and (3) the majority of the surface is covered 
in quiet chromosphere and a small fraction of moderately 
active regions. In the first case, as an active region on 
the star rotates into view, the equivalent width of Ca II 
K emission increases, while the Ha absorption line fills 
in and also increases in equivalent width (i.e. becomes a 
shallower and thus "less negative" absorption line). Once 
this isolated active region rotates out of view, both lines 
decrease. Our first case would show a positive correla- 
tion between Ca II K and Ha equivalent widths measured 
over time. In our second case, the majority of the star 
is producing both Ca II K and Ha in emission. As a 
region of quiet chromosphere rotates into view, Ca II K 
decreases slightly as both the basal and active Ca II K 
lines are in emission. The observed Ha line is now a su- 
perposition of an absorption line with the chromospheric 
emission line, causing the equivalent width to decrease. 
Again, one would observe a positive correlation between 
Ca II K and Ha. In our last scenario, the majority of 
the star is covered by basal chromosphere producing a 
weak Ha absorption line and Ca II K emission, and a 
moderately active region comes into view. If this moder- 
ately active region produces a deeper Ha absorption line, 
the equivalent width of Ha will decrease, while the Ca 
II K emission will be enhanced. In this case, one would 
observe a negative correlation between Ca II K and Ha. 
In between these extremes, variations in active region 
coverage and distribution across the stellar surface may 
explain much of the observed behavior between these two 
chromospheric tracers. 

5. COMPARISON TO UV AND X-RAY DATA 

In order to investigate the relationship between the 
chromosphere and the corona, we compare our optical 
measurements to archival near-UV and X-ray data. Of 
the 81 stars observed, 30 have archival X-ray measure- 
ments, and a subset of 7 of these stars have near UV 
(Mg II) measurements. All UV measurements come from 
active stars (those with Ha in emission), while the X- 
ray subsample also includes some moderate and weak 
activity stars (which have Ha in absorption). In Ta- 
ble 5, we provi de previously published Mg II and X-ray 
measurements (IWalkowicz et alj|2008b lBookbinderl ll985l: 
UohnsonlH98l IMathioudakis fc Dovleill989t 11992ft where 
available. 

Figure [T^] plots the X-ray flux versus the Mg II flux 
for this subsample. Both the X-ray and UV observa- 
tions were drawn from various sources: X-ray values were 
taken either from ROSAT observations of nearby stars 
(Huensch 1999) or EINSTEIN/EXOSAT values origi- 



nally reported in [Mathioud akis fe Dovld (|1989f ). noted 
in the plot legends as "MD89". Measurements of 
the Mg II flux were taken from either IUE observa- 
tions (IMathioudakis fc D oyle 1983) or ACS observations 
(jWarkowicz et al.ll2008f) . As some of the stars in this sub- 
sample have been observed by multiple instruments, the 
points in each plot are color coded by where the obser- 
vations originate, and multiple observations of a single 
object are connected by solid lines between points. Al- 
though a number of M dwarfs have been observed with 
newer X-ray instruments (such as XMM and Chandra) 
and with higher resolution in the UV (as with STIS or 
GHRS), our sample is co mprised only of M3 dwa r fs. We 
include the ROSAT and IMathioudakis fc Dovld (|1989l ) 
data sets because they are available for a number of stars 
in our sample, not just a select few of the most active 
ones (e.g., AD Leo and EV Lac). The solid line located 
to the right is the approximate power law fit to the Mg 
II and soft X-ray em ission in G and K dwarfs reported 
in lAvres et al.l (|1981ft . while the dashed line represents 
equality between the Mg II and X-ray emission. Our 
sample seems to obey a s imilar power law relation as 
that found by lAvres et alj (|1981ft for earlier type dwarfs. 

Our X-ray and Mg II detected subsample can also be 
compared to the Ha and Ca II K measurements reported 
here, as shown in Figures [T3] and Q3J respectively. In 
the plots comparing the flux in the Ha line, filled points 
indicate Ha in emission and unfilled points indicate Ha in 
absorption. The Ca II K line always appears in emission. 

Figure rjj] shows archival UV measurements of Mg II 
versus our measurements of Ha and Ca II K. As the Mg 
II observations are all of relatively active M dwarfs, this 
figure shows a positive correlation between Mg II and 
both Ha and Ca II K, reminiscent of the relationship 
between Ha and Ca II K in active stars shown in Figures 
4a and 5a. 

Figure [TBI shows a clear correlation between the optical 
chromospheric emission and the coronal X-ray emission. 
The relationship between the coronal emission and the 
Ha emission has somewhat less scatter than that between 
the X-ray and Ca II K flux (the three points that sit to 
the left of the main locus in Figure [T3a . which appear to 
have anomalously low Ha flux, are probably intermedi- 
ate activity stars whose Ha lines are almost completely 
filled in with chromospheric emission). The relationship 
between Ha and X-rays may be indicative of the link be- 
tween emission from the upper chromosphere (where Ha 
is formed) and the coronal emission. Figures [T2] and [Li] 
show that the Ca II K, Ha and Mg II chromospheric lines 
trace the coronal X-ray emission in the same way as one 
another, likely due to this subsample being comprised of 
active stars, whose greater activity leads to more emis- 
sion in all lines. 

Plasma temperature distributions have been calculated 
using Chandra data for only two of the most active stars 
in our sample, AD Leo and EV Lac; the bulk of the 
plasma lies between 2— 20 x 10 6 K, peaking arou nd 7—8 
x 10 6 K for both stars (|Robrade fc Schmittll2005T) . How- 
ever, ROSAT hardness ratios, which compare the flux 
in the 'soft' (0.14 - 0.42 keV; 'S') and 'hard' (0.52 - 
2.01 keV; 'H') bands, are available for a larger number 
of stars. The hardness ratio, (H— S)/(H+S), can be used 
as a rough proxy for coronal temperature, where harder 
emission (larger hardness ratio) corresponds to hotter 



coronal plasma. Figure [15] shows the hardness ratio for 
ROSAT detected stars versus the flux in Ca If K and Ha. 
It is evident from both of these plots that the chromo- 
sphcric emission, as traced by Ha is positively correlated 
with the hardness of the coronal emission, albeit with 
large scatter. In the case of Figure fT5b. the two points 
that appear to have low values of Ha flux for their coro- 
nal hardness are again the intermediate activity stars, 
where the equivalent width of the line is very small due 
to its being mostly filled in with emission. The Ca II K 
chromospheric emission also appears to be loosely corre- 
lated with the hardness of the coronal emission, but with 
far greater ambiguity than seen in the case of Ha. 

We can further examine the relationship between the 
coronal and chromospheric emission by looking at the ra- 
tio of the luminosity in Ha to the X-ray luminosity versus 
the hardness ratio. In Figure 1161 unfilled symbols indi- 
cate Ha in absorption, while filled symbols indicate Ha in 
emission. The two points that lie below this locus are the 
same two stars with very small Ha fluxes due to the line 
being completely filled by chromospheric emission. The 
median value for the sample as a whole is log(Lfr Q /La;) 
~ —0.6, which is cons i st wit h the values reported by 
iReid. Hawlev fc Mated (|1995l ) for low mass stars in the 
Hyades {\og{\j Ha /^x) ~ -0.6) and Pleides (log(Lff a /L x ) 
~ —0.8). For the most part, there is no obvious trend 
in the data- however for the most active stars, there ap- 
pears to be a slight indication that the X-ray luminosity 
increases proportionately with the coronal temperature, 
shown here as a slight negative correlation in the filled 
data points. Considering the dM and dMe stars sepa- 
rately, the median value of log(Lfj Q /La;) for the dM stars 
is log(Ljj a /L x ) ~ —0.45, somewhat higher than the dMe 
sample, where log^^/L^) ~ —0.84. These results are 
also c onsis tent with those reported in ISkumanich et all 
(|1984f) and|FJeminj (|1988f ) for field M dwarfs. 

Figure [TBI also compares the ratio of the chromospheric 
to coronal emission, respectively measured by Ha and X- 
rays, to model predictions of the minimum X-ray flux in 
late type stars. The colored bars in this plot correspond 
to our measured Ha flux divided by the range of X- 
ray flux generated by various models of coronal heating: 
model predictions of acoustically heated cor onae (log F x 
= 5.1 ergs s" 1 cm" 2 : iMullan fc 'Chendfl99l . and the X- 
ray flux predicted by models of magnetic field generation 
in the convective env elopes of non-rotat ing stars (log F^ 
= 3.8 ergs s _1 cm" - 2 : lBercik et~aLll2005[ ). Red bars indi- 
cate stars with Ha in emission, while blue bars indicate 
stars with Ha in absorption. Points lying below their ac- 
companying bars exceed the amount of coronal emission 
that can be explained by either of these models, while 
the X-ray flux of points lying within their colored bars 
could potentially be explained by either prediction. All 
of the active stars (dMe; filled diamonds) exceed the level 
of coronal emission produced in purely acoustic models, 
in accordance with the role global, rotationally depen- 
dent magnetic activity is believed to play in the heat- 
ing of the outer atmospheres of these stars. In the case 
of the dM stars (open diamonds), the behavior is more 
complicated- in most cases the coronal flux exceeds both 
predictions, but in some cases it does not. This obser- 
vation may provide another indication that two compo- 
nents are at work in the outer atmospheres of M dwarfs- 
hot, active regions associated with global magnetic activ- 



ity, coexisting with cooler basal regions that may be asso- 
ciated with small scale convectively-generated magnetic 
fields, analogous to quiet regions on the Sun. The rela- 
tive filling factors of these two components may explain 
much of the variation we see in M dwarf chromospheric 
and coronal emission. 

6. CONCLUSIONS AND FUTURE WORK 

The observations presented in this paper provide a rich 
set of empirical constraints for the construction of new 
model chromospheres. We find that for weak to inter- 
mediate activity stars, there is a range in the equivalent 
width of Ca II K emission for a small range of Ha ab- 
sorption. Our high resolution data indicate that the Ha 
line may undergo an initial increase in absorption, as 
Ca II K increases, with the line filling in and eventu- 
ally goin g into emission for the most active stars as ob- 
served bylStauffer fc Hartmannl (|1986h and predicted by 
ICram fc Giampapal ( 19871 ). From Ha observations alone, 
it is impossible to distinguish between a moderately ac- 
tive and an inactive star. Therefore, studies of activity in 
low mass stars that focus solely on stars with Ha in emis- 
sion are biased towards the most active chromospheres. 
These active stars have both Ca II and Balmer lines in 
emission, and more active stars have stronger emission 
in both lines. This positive correlation between Ca II K 
and Balmer emission in single observations of active stars 
does not imply, however, that the chromospheric tracers 
are positively correlated over time for every star. The 
variation in chromospheric emission of individual stars 
is most likely a function of the quiet and active region 
coverage unique to their particular stellar atmosphere. 

That lateral inhomogeneity may play a major role 
in the observed spectra of our M3 sample is not sur- 
prising, given observations of active regions on the 
Sun, and the inferred large filling factors of active re- 
gions deduced from M dwarf magne tic field measure- 
ments and observations of flares fe.g. iBasri et al.lll99Ct 
Uohns-Krull fc Valentil[l99l lO'Neal et al.ll2004h . In the 
scenario we propose, stronger chromospheric emission 
would correspond to a greater presence of active regions 
on the stellar surface due to the increased influence of 
the magnetic field. If Ca II K and Ha are produced pref- 
erentially by cool and hot components, respectively, one 
would expect that while greater chromospheric heating 
would cause more emission in both lines, the relation be- 
tween the two lines would be increasingly decoupled as 
the hotter active component becomes more prominent. 

We are in the process of further investigating whether 
two-component model atmospheres can account for the 
relationships between the Ca II and the Balmer lines. 
Using the data described here as empirical constraints, 
we are developing a new generation of quiescent non-LTE 
model atmospheres for M dwarfs. These two-component 
static models will be calculated using the n on-LTE ra- 
diativ e transfer program 'RH' described in lUitenbroekl 
(|2001h . RH allows for bound-bound and bound- free ra- 
diative transitions to overlap in wavelength and includes 
the effects of partial redistribution for strong bound- 
bound transitions, essential for proper treatment of the 
Ca II and Balmer lines. This effort complements recent 
work on modeling the outer atmospheres of M dwarfs 
during flaring states with the radiative-hydrodynamic 
non-LTE transfer code RADYN (Carlsson & Stein 1997; 
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Abbett & Hawley 1999; Allred et al. 2005). These new 
models will address both the observed line fluxes and cor- 
relations between Ca II K and the Balmcr lines obtained 
from our lower resolution data, and in addition will al- 
low us to compare the detailed model line profiles to our 
echelle observations for a subset of our sample stars. 
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Fig. 1. — Example HIRES (left) and DIS (right) observations of Ca II K and Ha for three stars of varying activity level. Top: AD Leo, 
high activity; Middle: Gl 362, moderate activity; Bottom: GJ 1125, low activity. It is evident that even stars that appear to be inactive 
at low resolution (e.g. GJ 1125) may still show chromosphcric Ca II K emission in high resolution spectra. 
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Fig. 2. — Comparison between Ca II K equivalent widths measured from low resolution data corrected for the photospheric contribution, 
and equivalent widths measured by fitting the Ca II K emission core in high resolution data. Diamonds denote measurements of our 
ARCES echelle data, while triangles indicate measurements from HIRES cchcllc data. The dotted line represents the locus of agreement 
for the two equivalent width measurements. For small equivalent width, the two methods of measurement produce comparable results, and 
the points lie along the line of agreement. However, for the most active the equivalent widths measured from the low resolution data are 
larger than those measured from the high resolution data. This difference is due to how the photospheric absorption line is accounted for 
in our two measurement methods: the method used for measuring the high resolution data underestimates the chromosphcric contribution 
to the line. The two points in (b) which lie significantly above the locus are two stars that were observed with DIS during higher activity 
states than their echelle observations. 
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Fig. 3. — Comparison between Ha equivalent widths calculated from DIS and echelle spectra. The dotted line represents the locus of 
agreement for the two equivalent width measurements. Diamonds denote measurements of our ARCES echelle data, while triangles indicate 
measurements of HIRES echelle data. As in the case of the Ca II K line, the two methods of measurement produce comparable results for 
small equivalent widths. However, for the most active the equivalent widths measured from the low resolution data are larger than those 
measured from the high resolution data. This discrepancy is caused by the Ha line falling close to the edge of the order in our echelle 
spectra- the continua in our echelle spectra have a different slope than the low resolution spectra, which causes the continuum value, and 
thus the equivalent width, to be underestimated in the echelle measurements. The point in (b) which lies significantly above the locus is a 
star observed with DIS during a higher activity state than its echelle observation. 




Fig. 4. — The total flux in Ca II K (red diagonal hatching) is comprised of a chromospheric emission core (blue horizontal hatching) 
superimposed over a photospheric absorption line (green vertical hatching). However, some flux in the inner line wings is also chromospheric, 
i.e. the flux exceeds that predicted by the radiative equilibrium model absorption line (green vertical hatching). Our method of measuring 
the Ca II K line corrects for the photospheric contribution to the line by subtracting off the radiative equilibrium profile, thereby measuring 
only the chromospheric flux in Ca II K. 



Ha Equivalent Width v. Call K Equivalent Width 

8[ 




Hex Equivalent Width v. Call K Equivalent Width 

0.4 L J 



0.2 

0.0 

•<. -0.2 

* -0.4 

-0.6 

-0.8 
-1.0 




Inactive 



Intermediate 
Chromospheres 



-1.0 -0.5 0.0 0.5 1.0 

EW C .„ K (A ) 

(b) 



1.5 



2.0 



Fig. 5. — (a): Total Ha equivalent width versus chromospheric Ca II K for the entire sample observed with DIS. For low to intermediate 
activity objects (those with Ha absorption) there is no obvious correlation with Ca II K. However, active stars show a positive correlation 
between the lines, (b): Total Ha equivalent width versus chromospheric Ca II K, shown at an expanded scale for the weak to intermediate 
activity stars. The perpendicular dotted lines plotted over the data indicate zero Ha and Ca II K equivalent widths- stars with Ha 
absorption (negative equivalent width) would traditionally be classified as "inactive" , but many of these stars possess small chromospheric 
Ca II K emission. Increasing absorption in Ha as Ca II K increases is not evident from these low resolution data. 
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Fig. 6. — (a:) Ha and Ca II K equivalent widths measured from our echelle spectra. Diamonds denote measurements of our ARCES 
echelle data, while triangles indicate measurements from our HIRES data. As in the case of the low resolution observations, Ha absorption 
is associated with a range in Ca II K emission, while Ha emission seems to be positively correlated with Ca II K. The weak to intermediate 
activity objects , plotted at an expand e d scal e in (b), seem to indicate an initial increase in Ha absorp tion with increasing Ca II K, as per the 
observations of Stauffcr & Hartmann (1986). The arrows indicate the schematic locus (predicted by Cram & Giampapa 1987) of an initial 
increase in Ha absorption as Ca II K increases, and subsequent filling in with emission as Ca II K increases even further. Interestingly, the 
locus appears to delineate the envelope of the observations. 
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Fig. 7. — Ha equivalent width versus the Ca II K surface flux measured from our low resolution DIS data. More active stars (those with 
larger Ca II K flux) have a wide range of scatter in their Ha equivalent widths. 





Fig. 8. — Equivalent widths for higher order Balmer lines versus Ca II K equivalent widths from DIS data. The higher order line and Ca 
II K show a similar relationship as Ca II K did with Ha in Figure 4(a). 




Fig. 9. — Multiple measurements of equivalent widths of the Balmcr lines and Ca II K in the active star AD Leo, measured from HIRES 
echelle data. Although a positive correlation between the instantaneous measurements of Ha and Ca II K is evident in Figures 2 and 3, it 
is clear that the two lines are not neccessarily positively correlated in time-resolved observations of a single star. 
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Fig. 10. — Continua near H« and Ca II K from our flux-calibrated DIS data versus absolute V magnitude. There is no disccrnablc trend 
between the continuum level and the absolute magnitude of the stars in our sample. 
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Fig. 11. — Equivalent widths of Ha and Ca II K (measured from DIS data) versus absolute V magnitude. The activity levels we determine 
for each star show no bias towards higher luminosity stars. 




Fig. 12. — X-ray flux versus Mg II flux for a subset of our sample. Different plotting symbols are used to indicate measurements by 
different instruments. Points corresponding to multiple observations of a given star by different instruments are connect ed by solid lines. Th e 
solid line at right represents the approximate power law fit to Mg II and soft X-ray emission of G and K dwarfs given in Avres et al.l (1198 lh , 
while the dashed line indicates the line of equality between the Mg II and X-ray emission. 
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Fig. 13. — X-ray flux versus the optical chromospheric flux as measured by Ha and Ca II K (measured from flux-calibrated DIS data). 
In (a), open circles correspond to Ha absorption, while filled circles denote Ha emission. As the Ca II K line flux we report is always due 
to the chromospheric emission core, all symbols in (b) are filled to indicate emission. The three points that lie to the left of the main locus 
in (a) are Ha lines with very low equivalent width, probably due to their being almost completely filled in by chromospheric emission. 
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Fig. 14. — UV chromospheric (Mg II) flux versus the optical chromospheric flux as measured by Ha and Ca II K. Different plotting 
symbols indicate Mg II measurements taken with either IUE or the ACS near-UV prism; multiple observations of stars are connected with 
solid lines. In (a), two of the stars with Mg II measurements have Ha absorption— these are indicated by open plotting symbols, while 
emission is denoted by filled symbols. In the Mg II sample of active stars, the positive correlation is reminiscent of that for Ca II K and 
Ha in our optical measurements. 
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(a) Hardness ratio versus Ha 



(b) Hardness ratio versus Ca II K 



Fig. 15. — Hardness ratio versus measurements of Ha and Ca II K flux from DIS data. Open diamonds in (a) denote Ha absorption, 
while filled diamonds indicate emission. Increasing hardness is indicative of hotter coronal temperatures, and correlates with increasing 
flux in the chromospheric lines. 
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Fig. 16. — Ratio of the luminosity in Ha to the X-ray luminosity versus the hardness ratio. Open diamonds indicate Ho absorption, 
while filled symbols indicate He* emission. Stars lie at a roughly constant value of hfja/^x ~ -0.6, consistent with observations of both 
field and cluster M dwarfs. 
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TABLE 2 

Wavelength Regions Used in Line Calculations 



Line Blue Continuum (A) Line (A) Red Continuum (A) 



H8 3850.00 - 3880.00 3884.15 - 3898.15 3910.00 - 3930.00 

Ca II K 3952.80 - 3956.00 variable 1 3974.80 - 3976.00 

H <5 4060.00 - 4080.00 4097.00 - 4110.00 4120.00 - 4140.00 

H 7 4270.00 - 4320.00 4331.69 - 4350.00 4360.00 - 4410.00 

H a 6500.00 - 6550.00 6557.61 - 6571.61 6575.00 - 6625.00 



"Integration region for the Ca II K line was chosen interactively for individ- 
ual spectra. 



TABLE 3 

Equivalent Widths Measured from DIS Data 



[RHG95] 


Name 


Ca II K 




He 


"He 


H8 
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<?H5 


H 7 




Ha 






(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 



460 GJ 3168 0.53 0.20 -5.00 0.18 -0.37 0.43 -1.00 0.27 2.24 0.12 -0.14 0.08 

492 Gl 109 0.39 0.21 -5.06 0.19 -0.53 0.45 -1.13 0.30 1.79 0.14 0.00 0.07 

836 Gl 179 1.06 0.29 -4.38 0.18 -0.56 0.41 -1.16 0.29 1.86 0.12 -0.23 0.11 

876 GJ 3333 0.46 0.15 -5.01 0.15 -1.48 0.38 -1.21 0.28 2.61 0.11 0.01 0.07 

919 Gl 203 0.12 0.29 -5.30 0.24 -0.72 0.41 -1.84 0.34 1.98 0.16 0.09 0.08 

933 Gl 206 8.55 0.36 7.55 0.21 4.12 0.36 2.97 0.33 8.78 0.14 4.57 0.10 

941 GJ 3356 0.50 0.30 -4.46 0.16 -0.68 0.41 -1.60 0.31 1.79 0.14 -0.09 0.08 

1068 GJ 3404A 0.18 0.29 -5.09 0.21 -0.54 0.45 -1.32 0.32 1.80 0.14 -0.08 0.08 

1089 GJ 3412 0.28 0.23 -5.46 0.20 -0.34 0.48 -0.94 0.29 1.62 0.13 -0.53 0.10 

1094 Gl 251 0.19 0.24 -5.01 0.19 -0.33 0.44 -1.28 0.29 1.97 0.13 -0.11 0.08 

1126 Gl 263 0.35 0.29 -5.05 0.22 -1.10 0.47 -0.30 0.32 2.09 0.15 -0.07 0.08 

1168 Gl 273 0.37 0.21 -5.27 0.17 -1.40 0.38 -1.11 0.30 2.38 0.12 -0.09 0.11 

1184 Gl 277B 5.11 0.16 -0.38 0.19 0.16 0.36 0.14 0.30 3.61 0.16 0.94 0.09 

1247 GJ 1105 0.62 0.23 -4.24 0.17 -1.55 0.39 -0.69 0.29 2.50 0.13 -0.07 0.09 

1328 GJ 2069A 16.91 0.29 16.82 0.18 8.77 0.25 7.14 0.27 13.02 0.13 6.02 0.10 

1353 GJ 3511B 0.26 0.30 -5.68 0.28 -1.06 0.56 -0.11 0.30 2.57 0.16 -0.12 0.06 

1405 GJ 3522 7.10 0.21 4.99 0.17 1.90 0.30 1.53 0.33 7.08 0.11 2.76 0.10 

1420 GJ 3527 0.15 0.25 -5.91 0.17 -1.78 0.38 -1.64 0.33 2.22 0.16 -0.22 0.09 

1441 GJ 3537 7.28 0.15 5.73 0.30 3.70 0.53 3.88 0.28 7.92 0.21 3.79 0.08 

1501 GJ 1125 0.25 0.28 -4.79 0.17 -1.31 0.40 -1.61 0.32 2.14 0.14 0.00 0.08 

1507 Gl 352A 0.19 0.19 -5.87 0.22 -1.37 0.40 -0.39 0.32 2.00 0.17 0.00 0.07 

1529 Gl 362 3.73 0.29 1.75 0.20 3.43 0.30 3.05 0.26 6.19 0.14 1.10 0.10 

1530 Gl 363 0.52 0.21 -4.28 0.19 -0.92 0.39 -0.48 0.27 2.40 0.13 -0.19 0.10 
1611 Gl 386 0.26 0.33 -6.07 0.24 -0.88 0.43 -0.09 0.32 2.25 0.18 -0.10 0.07 
1616 Gl 388 9.64 0.27 6.54 0.22 3.20 0.33 3.02 0.27 6.82 0.16 3.73 0.08 

1657 GJ 3612 0.03 0.27 -6.20 0.24 -0.88 0.44 -0.51 0.31 2.38 0.16 0.13 0.09 

1658 Gl 398 10.40 0.25 9.03 0.23 5.79 0.34 6.10 0.31 10.66 0.17 5.14 0.09 
1660 GJ 3613 9.86 0.25 7.09 0.25 3.76 0.38 3.41 0.27 7.38 0.19 3.96 0.09 
1690 G1 403 2.85 0.36 -5.15 0.24 -1.27 0.46 -0.49 0.33 2.45 0.16 -0.20 0.08 
1742 GJ 3647 9.41 0.32 9.25 0.27 7.10 0.38 6.43 0.36 11.08 0.19 5.39 0.08 
1795 GJ 3666 0.21 0.24 -5.72 0.24 -1.29 0.40 -0.39 0.32 2.24 0.11 0.05 0.09 
1838 G1 443 0.92 0.15 -4.89 0.24 0.07 0.45 -0.39 0.25 2.23 0.17 -0.22 0.08 
1845 G1 445 0.05 0.31 -6.10 0.24 -1.27 0.42 -0.46 0.31 2.25 0.13 -0.73 0.13 
1931 GJ 3719 9.95 0.21 8.29 0.25 6.88 0.47 5.37 0.26 10.16 0.16 4.83 0.09 
1988 G1 480 0.59 0.20 -5.28 0.22 -1.37 0.37 -0.34 0.32 2.05 0.16 0.05 0.08 
2009 G1 486 0.17 0.24 -5.68 0.32 -1.68 0.58 -0.36 0.33 1.75 0.16 0.01 0.09 
2088 GJ 3775 0.61 0.31 -4.13 0.17 -1.53 0.35 -0.12 0.31 2.29 0.15 -0.13 0.10 
2094 Gl 507B 0.26 0.33 -4.41 0.21 -1.56 0.39 -0.51 0.31 2.40 0.18 -0.13 0.09 
2175 GJ 3804 0.31 0.38 -4.42 0.21 -1.19 0.40 -0.34 0.30 2.57 0.20 -0.12 0.10 
2300 Gl 553.1 0.14 0.28 -4.79 0.19 -1.27 0.43 -0.56 0.29 2.53 0.15 -0.20 0.10 
2310 Gl 555 0.77 0.26 -3.41 0.17 -1.42 0.38 -0.39 0.28 2.31 0.14 -0.17 0.11 
2399 GJ 3892 0.40 0.29 -4.93 0.20 -1.37 0.34 -0.43 0.30 2.41 0.18 -0.24 0.09 
2410 Gl 581 0.08 0.32 -4.82 0.17 -1.58 0.38 -0.50 0.31 2.36 0.18 -0.17 0.08 
2456 GJ 3913 0.17 0.32 -5.03 0.19 -1.73 0.36 -0.54 0.30 2.66 0.15 -0.45 0.12 
2472 Gl 597 0.20 0.37 -4.81 0.18 -1.56 0.42 -0.74 0.32 2.43 0.25 -0.13 0.09 
2567 Gl 617B 0.63 0.28 -4.58 0.19 -1.35 0.36 -0.53 0.31 2.31 0.19 -0.27 0.08 
2587 GJ 3953 1.16 0.23 -4.00 0.18 0.06 0.40 -0.52 0.31 2.48 0.15 0.10 0.09 
2599 Gl 628 0.71 0.28 -4.86 0.17 -0.77 0.41 -0.42 0.28 2.20 0.10 -0.27 0.10 

2660 Gl 643 3.79 0.26 0.20 0.20 2.48 0.33 2.11 0.28 5.58 0.12 2.21 0.07 

2661 Gl 644A 0.44 0.51 -4.34 0.26 -0.58 0.45 -0.69 0.37 2.33 0.14 -0.14 0.10 
2665 GJ 1207 7.49 0.39 6.11 0.33 4.23 0.49 2.61 0.35 8.05 0.19 3.94 0.10 
2700 Gl 655 0.45 0.18 -5.10 0.18 -0.71 0.39 -0.45 0.28 2.11 0.13 -0.37 0.08 
2712 Gl 660A 0.56 0.24 -5.22 0.22 -1.29 0.44 -0.78 0.36 2.22 0.16 0.08 0.09 
2714 GJ 3992 0.47 0.17 -5.06 0.18 -0.92 0.37 -0.45 0.27 2.04 0.12 -0.41 0.10 
2717 Gl 661A 0.06 0.28 -5.88 0.21 -0.97 0.41 -0.60 0.29 2.50 0.16 -0.34 0.09 
2719 GJ 1212 0.22 0.32 -4.13 0.18 -1.55 0.39 -0.31 0.29 2.45 0.16 -0.20 0.09 
2744 Gl 669A 4.45 0.25 0.58 0.17 1.36 0.30 0.94 0.31 4.42 0.14 1.48 0.12 
2797 Gl 687 0.14 0.14 -4.77 0.27 -1.35 0.39 -0.76 0.41 3.00 0.15 -0.23 0.11 
2846 GJ 4038 0.42 0.18 -5.33 0.18 -0.67 0.41 -0.55 0.28 1.97 0.12 -0.16 0.08 
2849 Gl 699 0.23 0.46 -4.42 0.23 -1.52 0.44 -0.86 0.34 2.93 0.13 -0.78 0.13 
2851 GJ 4040 1.13 0.24 -4.19 0.20 -0.17 0.43 -0.16 0.33 2.30 0.15 0.15 0.08 
2888 GJ 4048 0.12 0.22 -5.88 0.30 -1.45 0.53 -1.07 0.34 1.56 0.18 -0.21 0.09 
2895 GJ 1226 0.13 0.20 -6.07 0.25 -0.68 0.47 -1.14 0.37 1.76 0.16 -0.09 0.09 
2898 Gl 712 0.09 0.34 -5.53 0.22 -1.34 0.38 -0.69 0.31 2.76 0.13 -0.19 0.10 
2901 GJ 4055 0.15 0.23 -5.42 0.23 -0.89 0.45 -0.79 0.33 1.97 0.15 -0.09 0.10 
2916 GJ 4062 0.66 0.22 -4.72 0.22 -1.22 0.46 -0.45 0.34 1.77 0.14 -0.13 0.09 
2921 GJ 4063 1.10 0.28 -4.24 0.18 -0.96 0.39 -0.33 0.28 2.20 0.13 0.08 0.09 
2924 Gl 720B 0.29 0.26 -5.04 0.20 -1.28 0.40 -0.43 0.30 2.10 0.12 -0.13 0.10 
2937 GJ 4068 5.76 0.39 3.74 0.24 0.98 0.33 2.21 0.33 6.45 0.17 2.82 0.09 
2940 GJ 4070 0.21 0.19 -5.50 0.19 -0.84 0.40 -0.43 0.28 2.18 0.12 -0.14 0.07 

2945 Gl 725A -0.12 0.45 -5.35 0.32 -0.79 0.56 -0.87 0.43 2.80 0.29 -0.29 0.11 

2946 Gl 725B 0.10 0.26 -5.98 0.21 -0.28 0.47 -0.58 0.32 2.28 0.13 -0.12 0.07 
2964 GJ 4083 0.16 0.24 -6.34 0.47 -0.84 0.45 -0.59 0.29 2.25 0.14 -0.19 0.09 
2976 Gl 735 5.76 0.26 1.52 0.17 1.67 0.27 1.07 0.30 4.59 0.14 2.29 0.08 
3013 GJ 4098 0.22 0.31 -4.89 0.18 -1.31 0.40 -0.44 0.30 2.49 0.16 -0.49 0.10 
3032 GJ 9652A 10.77 0.27 9.09 0.22 4.26 0.30 4.19 0.29 9.71 0.18 5.40 0.11 
3058 GJ 4110 5.84 0.28 2.29 0.20 2.06 0.29 2.29 0.32 6.35 0.14 2.74 0.11 
3306 GJ 9721B 3.73 0.59 3.26 0.37 -0.17 1.05 2.33 0.39 6.38 0.31 2.69 0.09 
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Equivalent Widths Measured from Echelle Data 
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TABLE 5 

Stars with UV and X-ray Measurements 
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